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Abstract. We map the problem of the orbital excitation (orbiton) in a 2D antiferromagnetic 
and ferroorbital ground state onto a problem of a hole in 2D antiferromagnet. The orbiton 
turns out to be coupled to magnons and can only be mobile on a strongly renormalized scale by 
dressing with magnetic excitations. We show that this leads to a dispersion relation reflecting 
the two-site unit cell of the antiferromagnetic background, in contrast to the predictions based 
on a mean-field approximation and linear orbital-wave theory. 



1. Introduction 

Calculating collective excitations in nonfrustrated but interacting spin systems may be 
considered as one of the few extremely successful stories in solving the strongly correlated 
electron systems. Independently on whether the approach which is used to calculate such 
excitations stems from the linear spin wave theory (for ordered states) or the Bethe Ansatz 
solutions (for quantum disordered states in ID), a remarkable qualitative and quite often also 
quantitative agreement with a number of experimental results is usually achieved, see e.g. [T]. 

On the other hand, collective excitations in interacting orbital systems have attracted much 
less attention of the community, primarily due to experimental problems in detecting the 
collective orbital excitations (orbitons) [2j. A priori, excitations in the orbital degree of freedom 
of Mott insulators should behave in a similar manner as those of the spin sector: in both cases 
the large Coulomb repulsion in the Mott insulators 'freezes' the charge degrees of freedom and 
largely reduces the effective low energy Hilbert space so that instead of electrons it is given by 
spin or orbital (pseudospin) degrees of freedom. In place of the 'spin-wave theory' for magnetic 
excitations, one would then use 'orbital- wave theory' [3]. 

It has recently been shown [3], however, that this orbital- wave picture does not hold in 
many realistic systems. Here, we prove that in a generic 2D Kugel-Khomskii (KK) spin-orbital 
model [5] with the antiferromagnetic (AF) and ferroorbital (FO) ground state the linear orbital 
wave theory breaks down and one needs to explicitly take care of the coupling between spins 
and orbitals. Furthermore, we show that for realistic parameters of the KK Hamiltonian, the 
dispersion of the coherent excitation changes qualitatively and reflects the two-site unit cell 
of the antiferromagnetic background, in contrast to expectations based on linear orbital-wave 
theory. 
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Figure 1. Spectral function for 
the orbital excitation 0(k, oS) [Eq. 
([21)] calculated in the linear orbital 
wave picture [(a) and (c)] and using 
the exact mapping onto the effec- 
tive t-J model [(b) and (d)]. Upper 
(lower) panels for equivalent (in- 
equivalent) values of the hoppings 
between the two orbit als, see text. 
Note the different scales on upper 
and lower panels and that the in- 
coherent spectrum in panel (d) is 
extremely small and not visible at 
this scale. Parameters: E z = 5J 
and broadening r] = 0.05J. 



2. Model and results 

The Hamiltonian of the generic 2D KK model that we want to study describes a two-orbital 
system with orbital-conserving hopping and reads 
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where S (T) are the spin (pseudospin) operators which fulfill the SU(2) algebra and J is the 
superexchange constant which describes the interaction of S = 1/2 spins and T = 1/2 (orbital) 
pseudospins. Besides, E z is the crystal field which is implicitly assumed to be so large that in 
the ground state \<j)) of H only one type of orbital is occupied (the FO order of orbitals) and 
consequently there is long-range AF order in temperature T = 0. 

The purpose of this study is to calculate the orbital excitations from the ground state \(p) of 
Hamiltonian ([1]), as defined by the spectral function 0(k,uj) : 
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where E& is the ground state energy and is the momentum-dependent orbital raising 
operator. 

In the 'standard' orbital- wave approach [6, 7J, this spectral function can be easily calculated 
(see, e.g., Ref. [1]) by mean- field decoupling of spins and orbitals and by: (i) introducing the 
Holstein-Primakoff bosons for orbital pseudospins (cf. Ref. [8]) and (ii) neglecting interactions 
between bosons. We then obtain: 



0(k,w) 



u - E z + izJow(l - 7k) 



(3) 



with z = 4 is the coordination number in 2D and 7k = (cos k x + cos k y ) /2 is the lattice structure 
factor in 2D. Here J w = 2J(^[Sj • Sj + \\4>) ~ -0.16J where \cj)) is the 2D AF, i.e. spin-only 
part of ground state \4>). The spectrum consists of a single quasiparticle peak with a relatively 
large dispersion, cf. Fig. [Ha) f° r ^ ne spectrum along the high symmetry direction (0, 0) — > (tt, it) 
in the Brillouin zone. 
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Figure 2. Artist's view of the defects 
(springs) in the spin background (up 
and down arrows) created by the moving 
orbiton (bold oval). Only if these defects 
are healed by quantum fluctuations the 
orbiton can move coherently - otherwise, 
the orbiton is trapped in a string-like 
potential, see text. 



However, it was shown in detail in Ref. [1] that the above mean-field decoupling of spins 
and orbitals does not lead to a reliable description of orbital dynamics in a coupled spin-orbital 
system. Therefore, following Ref. [4j and using inter alia the Jordan- Wigner transformation [9] 
for spin and orbital degrees of freedom^ we map the above problem of a single orbiton in the 
AF and FO ground state Eq. ([2|) onto a a problem of a single hole in the AF ground state 
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described by an effective t—J Hamiltonian 
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Here the electron operators p- lcr (p\ a ) act in the restricted Hilbert space without double 
occupancies and create (annihilate) holes in the 2D AF ground state \4>). Besides, n\ = J2 a n 'w 
and the effective hopping parameter t describing the motion of orbital excitation is given by 
t = J/2. 

The spectrum Eq. was calculated using the self-consistent Born approximation on 32 x 32 
sites (which gives results in very good agreement with the exact diagonalization studies, cf. Ref. 
|10| ) and consists of a large quasiparticle peak and a very small incoherent spectrum, cf. Fig. 
[T^b). The bandwidth W of the coherent orbiton peak, as obtained from Fig. [T](b) or from 
perturbative calculations cf. Ref. [10], is W ~ 2t 2 / J ~ J/2 (for t = J/2) and thus somewhat 
smaller than the mean-field one for which W ~ 0.64J. In more striking contrast to the mean- 
field results, the dispersion has a minimum at (tt/2,tt/2) instead of (0,0), as a result of the 
coupling between orbitons and the AF 'background'. 

In many realistic cases, the spin-orbital model JT]) has to be modified to account for the 
different hoppings for electrons in the two active orbitals. We have verified that the effective 
hopping t is then no longer t = J/2 but instead t = 2t%t-2,/V where t\ and £2 denote the bare 
hopping of an electron in orbital 1 and 2 and U is the on-site Hubbard repulsion. Furthermore, 
quite often the hopping of the occupied orbital (defined below as orbital 1) is larger, as the 
gain in hybridization energy contributes to its energy gain with respect to the empty orbital, 
i.e., t\ > t%. For example for copper oxides with an occupied x 2 — y 2 orbital and empty tig 
orbital, ti ~ t\/2 and thus t = J/4 since J = At 2 /U. In the orbital -wave picture, the dispersion 
width is then likewise reduced by a factor of two to W = 0.32J, but the effect is stronger in the 
interacting t-J model, where we find W ~ 0.125J, cf. Fig. QJc)-(d). 



1 See also supplementary material in Ref. [3] which discusses the validity of this mapping in the 2D case. 



3. Discussion 

Let us now try to understand the physical origin of these large differences between the orbiton 
motion obtained in the mean-field approximation and in a far more exact way (i.e. via mapping 
onto the effective t—J model), cf. Fig. [2j As in the mean- field case the orbiton moves freely 
and does not excite any spin excitations, the whole effect of the spin background is taken into 
account by renormalizing the superexchange constant J into its reduced value Jow- This stays 
in contrast with the result obtained from the mapping onto the t—J model: now the orbiton is 
effectively a 'hole' in the AF spin background which, when moving, flips the spins of the 2D 
AF, cf. Fig. [21 In fact, the further the orbiton moves the more spin flips (defects in 2D AF) 
it excites which leads to the so-called string-like potential acting on the orbiton, cf. Fig. [21 
The orbiton can then move coherently only due to the presence of quantum spin fluctuations 
in the ground state which 'heal' these defects, 'attach' to the orbiton and jointly move as an 
orbiton quasiparticle on a renormalized scale (cf. Ref. |1R [TO] with the hole in 2D AF moving 
effectively as a spin polaron). We also note that this result is entirely different from that in 
the ID case [3J, where the orbital flip excitation fractionalizes into free 'spinon' and 'orbiton' 
quasiparticles, similar to decay of a doped hole into spinon and holon. 

4. Conclusions 

In summary, we showed that a single orbital excitation in the 2D AF and FO ground state moves 
similarly to a hole in the 2D AF state. Therefore, it can only be mobile by coupling to the spin 
fluctuations which leads to a minimum in the orbiton dispersion for the (it/2,tt/2) point in the 
Brillouin zone and to the very strong renormalization of the orbiton bandwidth. In particular, 
in the realistic case which can correspond to the dd excitations in copper oxides, the orbiton 
bandwidth is of the order of ~ 0.1 J which explains why typically such dd excitations in the 
AF ground state are nearly dispersionless and can be well modelled using quantum chemistry 
calculations on small clusters [T2"] . 

Acknowledgments 

We acknowledge support of the Alexander von Humboldt Foundation (K.W.) and the DFG 
Emmy Noether Program (M.D.). 

References 

[I] Coldea R et al. 2001 Phys. Rev. Lett. 86 5377 
Walters A C et al. 2009 Nat. Phys. 5 867 

[2] Gmninger M et al. 2002 Nature 418 39 

Ulrich C et al. 2009 Phys. Rev. Lett. 103 107205 
[3] van den Brink J, Stekelenburg W, Khomskii D I, Sawatzky G A and Kugel K I 1998 Phys. Rev. B 58 10276 

van den Brink J, Horsch P, Mack F and Oles A M 1999 Phys. Rev. B 59 6795 
[4] Wohlfeld K, Daghofer M, Nishimoto S, Khaliullin G and van den Brink J 2011 Phys. Rev. Lett. 107 147201 
[5] Kugel K I and Khomskii D I 1982 Sov. Phys. Usp. 25 231 
[6] Khaliullin G and Oudovenko V 1997 Phys. Rev. B 56 14243 
[7] Oles A M, Khaliullin G, Horsch P and Feiner L F 2005 Phys. Rev. B 72 214431 
[8] Wohlfeld K, Oles AM and Horsch P 2009 Phys. Rev. B 79 224433 
[9] Jordan P and Wigner E 1928 Z. Physik 47 631 
[10] Martinez G and Horsch P 1991 Phys. Rev. B 44 317 

[II] Schmitt-Rink S, Varma C M and Ruckenstein A E 1988 Phys. Rev. Lett. 60 2793 
[12] Hozoi L, Siurakshina L, Fulde P and van den Brink J 2011 Sci. Rep. 1 65 



